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Incorporating many Distributed Generators (DGs) technologies in power 
system networks has grown rapidly in recent years. Distributed generation 
(DG) plays a key role in reducing power loss and enhancing the voltage 
profile in radial distribution networks. However, inappropriate DGs site or 
size may cut network efficiency; moreover, injecting harmonics is one of the 
integration concerns of inverter-based DGs. Two-procedure based approach 
is introduced in this paper. The first procedure solves the DGs siting and 
sizing problem, as a multi-objective one by improving the voltage profile of 
the whole distribution network and also reducing its power loss. A weighted 
sum method is presented to create the Pareto optimal front in this procedure 
and get the compromised solution by applying a novel metaheuristic 
optimizer, named Crow Search Algorithm (CSA). A modification on CSA is 
also proposed and applied to improve its performance. The achieved solution 
for inverter-based DGs placement and size is checked in the second 


procedure to make sure the accepted voltage THD at all buses by 
implementing detailed simulation for the tested system using 
Matlab/Simulink. The proposed approach has been tested on IEEE 33-bus 
radial distribution system with photovoltaic DGs. To confirm the superiority 
of the modified CSA algorithm in terms of quality of solution, its achieved 
results are compared with the results offered by the original CSA algorithm 
and published results of some other nature-inspired algorithms. 
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1. INTRODUCTION 

Radial distribution systems (RDS) are commonly implemented in rural or suburban areas for their 
simplest structure and low initial installation cost as they are fed from utility grid at only one end. However, 
they suffer from several disadvantages. Their main disadvantages are the heavy loaded distributor at its end 
near to the substation, the customers at the distant end of the distributor face serious voltage fluctuation 
problem and in addition that distributor faults will cause interruption in supply to large number of consumers 
connected to that distributor. Actually, the overall power losses within the transmission and distribution 
systems limit the capability of the radial distribution system. In fact, distribution power loss accounts for the 
greatest amount of loss as 13% of the generated energy, which is dissipated as a heat within conducting 
materials [1]. This non-negligible quantity of losses results in a direct economic impact and also affects the 
overall capability of RDS. 

DGs at the distribution side of the network are projected to play a significant role in satisfying the 
rapidly increasing world’s electricity demand. DGs differ with respect to purpose, voltage level, size, 
technology and others [2]. In general, DG size ranges from some kW to hundreds of MW. It can be classified 
by micro DG for less than 5 kW; small DG from 5 kW to less than 5 MW; medium DG from 5 MW to less 
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than 50 MW and finally large distributed generation from 50 MW to less than 300 MW [3]. Based on 
technology and connection with distribution network, DG can also be classified as reciprocating engines, 
micro-turbines, combustion gas turbines, fuel cells, photovoltaic systems, small hydro turbines and wind 
turbines [4]. The main advantageous applications of DG can be summarized as less expensive standby 
generation instead of grid extension and improving system performance. Besides, DG units using renewable 
energy resources offer environmentally friendly solutions than the traditional generation by avoiding the 
large measure of CO2 emissions. On the other hand, most DG technologies use power electronic converters 
to manage the ability of these generators. Such converters negatively affect the grid harmonics level that need 
to be carefully monitored to guarantee the allowable limits and ensure satisfactory operation of 
the system [5]. 

Determining the suitable place and size of DG is vital to get greatest benefits from integrating DG 
with distribution systems. Inappropriate sitting may scale back the advantages and drive to overall poor 
system performance [6]. Several approaches are proposed and developed in literature to decide the right site 
and size of DG units achieving one or several goals with or without network reconfiguration. The analytical 
and improved analytical methodologies for minimizing power loss have been introduced in [6-7]. A power 
stability index is employed in [8] to decide the more sensitive bus for DG unit placement, and an exploration 
rule was introduced for getting the optimum size to reduce power loss. Several researchers have also applied 
different metaheuristic techniques for getting the optimum location and size of DG units. In [9], the genetic 
algorithm (GA) is primarily used for solving multi-objective optimization for the most effective size and 
allocation of DG units with multi-system constraints whereas in [10-11], the particle swarm optimization 
(PSO) has been implemented to select the optimum size and placement of one DG unit. Both GA and PSO 
are joined in [12] to improve their performance for determining the DG optimum size and placement. 

In recent years, there has been an increasing amount of literature on applying very recent developed 
optimization tools for optimally DG allocation such as cuckoo search, bacterial foraging optimization 
algorithm, flower pollination algorithm, backtracking search optimization, ant colony optimization and 
teaching learning-based optimization [13-18]. During most of these optimization approaches, the generation, 
constraints of bus voltage and branch loading are incorporated in the fitness function as penalty factors to 
meet a valid solution such that all optimal power flow variables stay within their permissible limits. 

In the recent decade, the new Crow Search Algorithm (CSA) is applied in several research studies to 
get the optimal size and site of DGs [19-23]. The simultaneous optimal placement of DGs and capacitor units 
in RDS is studied in [20]. Reducing the power loss and improving the voltage profile of the RDS using 
multiple DGs were the goals of [21, 23], while maximizing the overall saving and reducing the network 
losses were the main objectives of [22]. Obtaining the optimal size and site of DG and DSTATCOM in 
practical distribution system, to reduce the distribution system loss, is fully introduced in [19]. 

In present work; modified CSA is implemented to formulate multi-objective function for getting 
optimal DG size and allocation. It is applied for single DG or multiple DGs units to reduce the active power 
losses and improve the voltage profile in RDS. The proposed technique will be tested on IEEE 33-bus test 
radial system. This paper will be organized as follows: problem formulation is presented in Section 2, 
while in Section 3 the metaheuristic CSA is introduced with the proposed changes. In Section 4, detailed 
simulation results are discussed. CSA algorithm and the modified one are compared against some other 
nature inspired algorithms in Section 5; and finally the conclusions are drawn in Section 6. 


2. OPTIMAL PLACEMENT AND SIZING PROBLEM FORMULATION 
This section will discuss briefly the applied method for load flow analysis while the mathematical 
formulation of the objective function for optimum DGs placement and sizing will be introduced in details. 


2.1. Load Flow Analysis 

The simple forward-backward sweep method introduced in [24] is applied here to solve the RDS 
load flow analysis and calculate the overall power loss and bus voltages at all buses. It begins at the end bus 
and precedes backwards the source node. The convergence is checked using the node voltage calculated in 
backward sweep. If the obtained voltage has less difference than the convergence criterion, the process will 
end, but if the voltage does not meet convergence limit, forward sweep begins to compute the voltages at 
each bus starting from the feeder source bus. 

According to the small sample network shown in Figure 1, the current in the branch section line I,, 
and the voltage V, at node ‘y’ are computed according to (1), known the load active and reactive powers 
(P;,- Qi) at node ‘x’ & the branch resistance and reactance (R, . X,.) respectively and the voltage V, at node 
“x’. Real power loss in the branch section connecting nodes x and y is expressed in (2): 
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Figure 1. Sampled distribution system for load flow analysis 


2.2. Mathematical Problem Formulation 
2.2.1 Multi-Objective Approach for Optimal DG Sizing and Placement 

As a matter of fact, the bus voltage of a typical RDS often experiences fluctuations and may 
collapse under some critical loading conditions with increased load demand. Therefore, the main purpose of 
the present work is to get the optimal location and size of multiple added DGs in RDS while achieving 
minimum real network loss and improved voltage profile. Accordingly, the objective function can be divided 
into two parts: power loss minimization and voltage deviation minimization. These two single objective 
functions are merged to form a multi-objective optimization problem, whose (OF) is the cost function. As 
given in (3), the overall multi-objective function is expressed using a linear combination of the two 
objectives OF, and OF, via the weighted sum method, where w,, w> are the weighting factors for OF, and 
OF, respectively. 


OF = w, x OF, + w, x OF, (3) 


Assuming that the studied RDS that has a rated network voltage V;ateq (equals 1.0 pu), consists of 
number of buses equals No_b and number of branch lines is No_lines. Each branch j connecting two buses 
has a branch resistance Rj, carries a current with the value I, and the voltage at the iy, bus is denoted by Vj. 
Therefore the two objectives OF, and OF, are expressed by (4) and (5). As clearly shown, the total active 
power loss of the system (Pr joss) can be calculated based on the summation of the power loss in all RDS 
branches, while the overall voltage deviation (VD-) is calculated based on the summation of the square value 
of bus voltage deviation (VD;) at all buses. 


OF, = minimizing total active power loss (Pr joss) = min Bee I? R; (4) 


OF, = minimizing overall voltage deviation (VD) = min N°? VD, = min YN?” (V, — Veatea)” 


(5) 


In fact, the weighted sum method is applied here to create the Pareto optimal front. This might be 
noted that each objective functions in (3) is normalized via dividing it by its base value which makes the 
objective function dimensionless and normalized and that prevents any scaling problem as discussed in [25]. 


2.2.2 Constraints 

The applied formulation as a multi-objective function is subjected several constraints: network 
operating constraints and DG constraints. The operating constraints define the voltage limits at all network 
buses, line capacity limits of different lines and the active power balance equations as follows: 


Vin < Vj < Vmax fori = 1to No_b (6) 


Ty} < Imax(j) for j = 1 to No_lines (7) 
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Psupply oH yes Poc k = Poaa + Pr joss (8) 


Where: 

Vin and Vinax are the smallest and the greatest allowable voltage V, at each bus i. Such limits are 
expressed according to [26] by + 5% of the rated network voltage V-ateq for all network buses. 

Tmax (Jj) is the maximum current carrying capacity [26-27], the current J; flowing in the line j should 
be less than Imax (j)- 

Psupply 18 the active power supplied by the grid as the original primary source, while Poag is the total 
active load of the system. Ppg , is the active power penetrated by k,, DG unit for N no. of DG units. 

On the other hand, some constraints about the size of each added DG at k,, bus (Ppg x) and the 
overall DG penetration should also be considered as follows: 


Pog min = Poc k s Pog max (9) 
Pr_pG min < diet Pock < Pr_pG_max (10) 
Where: 


Pog min 2Nd Ppg max denote the minimum and maximum allowed output active power of ky, DG 
unit respectively, 

Pr pg min@nd Py pg max are the total minimum and maximum allowed output of all added DG units 
(N units), assuming that all DG units generate active power only. 


2.2.3 Ensuring the voltage THDv at All Buses 
Distorted currents absorbed by harmonic-producing loads also distort the supply voltage as they 
pass through system impedance. Therefore, a distorted voltage can be presented to other end users. The 
power quality industry has developed certain indices to assess service quality related to distortions caused by 
harmonics. The total harmonic distortion is one of the most commonly used indices for measuring the overall 
waveform harmonic content taken into account the contribution of all harmonics components [28]. IEEE Std. 
519™- 2014 is developed for utilities and customers to limit the harmonic content, the boundaries on the 
foe) 
h=2 Vi 
total harmonic distortion on the output voltage (THD, =——-~—————) are listed in Table 1 for 
Fundamental voltage 
distribution networks levels. 
In the second procedure of the proposed approach, the impact of DGs on harmonics is investigated. 
The achieved former solution for DGs placement and size is tested to check THDy at all buses by simulating 
the RDS using Matlab/Simulink where the installed DGs type is assumed inverter-based. When the limits are 
still confirmed, the achieved solution is accepted. However, if the limits are violated, another solution is 
selected or inserting a mitigation filter is suggested. 


Table 1. Voltage Distortion Limits in Distribution Networks According to IEEE Std. 519™-2014 
Bus voltage V at point of common coupling _ Individual Harmonic (%) Total Harmonic Distortion (%) 
V<10kV 5.0 8.0 
LkV<V<69kV 3.0 5.0 











3. MODIFIED CROW SEARCH ALGORITHM (CSA) OPTIMIZER 

Crow search algorithm is a recent metaheuristic algorithm developed by Askarzadeh [29], inspired 
on the intelligence behavior conducted by crows of hiding their excess food in a place and get it back when 
needed. CSA has been applied for different problems with different constraints. As an algorithm based on 
population, the size of the flock is confirmed by N¢ individuals (crows) which are of n — dimensional where 
n denotes the problem dimensions. Each crow (individual) is assumed to have the capability of remember the 
best visited location to hide food. The position is then modified according to Pursuit & Evasion behaviors. 
Pursuit: a crow j follows crow i with the purpose to discover its hidden place. The crow i does not notice the 
presence of the other crow, as consequence the purpose of crow j is achieved. Evasion: the crow i knows 
about the presence of crow j and in order to protect its food, crow i intentionally takes a random trajectory. 
This behavior is simulated in CSA through the implementation of a random movement. The type of behavior 
considered by each crow i is determined by an awareness probability (AP). The flight length (FL) parameter 
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indicates the magnitude of movement from crow position towards the best position of crow j [30-31]. In fact, 

CSA consumes less time than other metaheuristic methods as it has fewer parameters to adjust: AP and FL, so 

it is easy to be implemented. To offer a good balance between diversification and intensification in CSA, the 

AP can control both factors. When AP decreases, the algorithm tends to search in local optima therefore, the 

intensification increases. When AP increases, the algorithm tends to search in global optima therefore, the 

diversification increases. 

In [19-23], the original CSA is applied to get the optimal size and site of DGs where the variable 
represents the position of the crows in the space is assumed randomly based on lower (Lg) and upper (Ug) 
bounds of the variable. A modified CSA is proposed in this article, where the CSA variable is calculated 
depending on Gaussian and Cauchy density functions [32], which are shown in Figure 2. 

a) The Gaussian function or the normal distribution is defined as a continuous function that approximates 
exact binomial distribution for the events. It provides probability for some real observations to fall 
between any two predefined real numbers or limits as the function tends to zero on either side. Its curve 
resembles the shape of a bell, thus it is informally known as the bell curve. For the mean or expectation of 
the distribution (1) and the variance (07), the probability density of the normal distribution is expressed 


1 Ge—w? ; ; Oey 
by F(x) = 5 € 207 . In paper work, the variance is assumed to be unity (i.e. o7=1). 
210 


V CTU 
b) The Cauchy distribution is a continuous probability distribution known as Cauchy—Lorentz distribution or 
Breit—Wigner distribution. It is expressed by: F(x) = Ser for t > 0 as a scale parameter. 
x 





— Gaussian 
-- Cauchy 


1m (1+(LB+rand*(UB-LB)-(LB+UB)/2) 2 


1 2% at) Brg (LB Hand" (UB-LE)(LB+UB)/2)/2 0? 




















LB (LB+UB)/2 UB 
Figure 2. Gaussian and Cauchy distribution functions 
In order to check the performance of the modified CSA on larger-scale optimization problems, four 
well-known benchmark functions, discussed in [33], are solved in 10 dimensions using CSA (original 


algorithm), CSA-G (based on Gaussian distribution), and CSA-C (based on Cauchy distribution). The 
obtained results are summarized in Table 2. As shown, CSA-G outperforms the other methods. 


Table 2. Performance of Modified CSA for Benchmark Test Functions 








Function Name Function Formulation Search Space Optimal Value CSA CSA-G CSA-C 
n 
Sphere Function y xX? [-100,100] 0 9.54 e-13 1.39 e-36 2.69 e-35 
i=1 
n 
Rosenbrock Function yy (100(xj41 — X7)? + (x; - D*) [-30,30] 0 1.52 0.75 0.77 
i=1 
1 n n xi 
Griewank function —) (X))? - I] cos(—) +1 [-600,600] 0 0.0099 0.0059 0.229 
4000 i=1 i=1 Vi 
n n 
Schwefel function y¥ |x;| - I] |x;| [-10,10] 0 9.37 e-6 3.29 e-18 1.13 e-17 
i=1 i=1 





4. TEST SYSTEM AND SIMULATION RESULTS 

The system under study is IEEE 33-bus radial distribution network shown in Figure 3. Without DG 
installation, it has a total load of 3.72 MW, 2.3 MVAR and real power loss of 211 kW. All system buses are 
not within the limits + 5% of the rated network voltage and the least voltage is 0.9037 pu at bus 18 [34]. 

In this paper, two-procedure based approach is introduced. The first one solves the DG siting and 
sizing problem as a multi-objective one by improving the voltage profile of RDS and reducing its power loss 
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as two conflicting functions using a modified CSA-G. Then the solution is checked in the second procedure 
to ensure accepted THD, at all buses. In the following sections, achieved results are introduced. The results 
for solving optimal DG placement and sizing problem using multi-objective optimization approach are 
introduced in first part. The impact of DG on THD, is investigated in the second one. 


23 24 25 






29 30 31 







32 33 


20 21 22 


Figure 3. IEEE 33-bus radial distribution network 


4.1. Results of Optimal DG Sizing and Placement 
Installation of single and multiple DGs are considered in the tests, which ended with a comparative 
study of all cases. 


4.1.1 Case 1: Installing Single DG 

In this case, a single DG unit is assumed to be integrated to the 33-bus RDS. By applying different 
values for weighting factors w,, Wz, sets of solutions are achieved using CSA-G. Figure 4-a) shows the 
solutions set or the Pareto front which consists of 10 trade-off solutions (values of w,, W2, change from 0.1 
to 0.9 with a step of 0.1). The point of intersection of the two curves denotes the compromised solution. 
The Pareto front for such optimal solution ensures that a DG size of 3114.5 kW can be installed at bus 7 to 
give the least voltage deviation of 0.02342 pu, and lowest reduction in real power loss from 211 kW without 
DG installation to 118.5 kW with a reduction of about 43.83% at weighting factors of 0.47 & 0.53. 


4.1.2 Case 2: Installing Multiple DG units 
In this case, the multi-objective optimization problem is resolved by installing two DG units 
simultaneously to the tested RDS and three simultaneous DG units. The Pareto front curves using the 
modified CSA-G are shown in Figure 4-b) and Figure 4-c). 
a. Point B represents the compromised solution for adding two DG units of size 1344.7 and 1161.3 kW at 
bus 12 and 30 respectively to provide the minimum voltage deviation at 0.007107 pu, and real power 
has reduced to 94.18 kW (55.36% reduction without DG installation) at weighting factors 0.52 & 0.48. 
b. Point C characterizes the compromised solution for adding three DG units of 846.2, 931.6 and 1296.2 
kW to be installed at buses 13, 24 and 31. Such sizes ensure the minimum voltage deviation at 0.0066 
pu, and reduction of about 62.5% in real power loss (79 kW). The weighting factors w,, wz are 0.58 & 
0.42 respectively in this case. 


From the quick comparison of the Figure 4 graphs, it can be concluded that multiple DG 
installations will reduce the real power losses more than using a single DG unit and also the voltage deviation 
is decreased significantly. Actually, the size of the added DGs regarding the reduced loss reduction can be 
economically evaluated by a simple economical index. Besides, the effect of adding DGs on the voltage of all 
buses can be assessed using the index of the average voltage deviation. 

For further evaluate the proposed changes on CSA, the three compromised solutions A, B and C, 
are fully summarized for all applied methods of CSA (original CSA, CSA-G, and CSA-C) in Table 3 
according to such proposed indices: 

N 

Economical Index = >— Biss “G.k (11) 

losswithout DG ~ T-losswith DG 


No 
Pao > ee by -<Vi wi 2 
Voltage Index (p. u) = Average deviation per bus = ——— (12) 
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(c). Three DG units 


Figure 4. Compromised solutions using modified CSA-G 


Table 3. Comparison of the Three Compromised Solutions (A, B And C) for Single and Multiple DG Units 














No. of Applied DG Size @ Bus ey Pr toss Ec. Voltage Lowest bus voltage 
DG units Algorithm (kW) pees (kW) Index Index (pu) (pu @ bus) 
CSA 2534.7 8 0.42, 0.58 132.1 32.3 0.022804 0.9532 @ bus 18 
One CSA-G 3114.5 7 0.47, 0.53 118.5 33.6 0.02664 0.9534 @ bus 18 
CSA-C 2579.3 6 0.61, 0.39 153.3 44.7 0.02009 0.9534@ bus 33 
CSA 1073.5, 1326.2 16,26 0.48, 0.52 103.9 22.4 0.021105 0.9537@ bus 33 
Two CSA-G 1344.7, 1161.3 12,30 0.52, 0.48 94.18 21.45 0.01468 0.9752@ bus 33 
CSA-C 2494.3, 811.7 28,16 0.35, 0.65 121 36.73 0.00798 0.982@ bus 25 
CSA 1251.3, 1271.9, 611.6 14,2430 0.38, 0.62 87 25.28 0.01834 0.9601 @ bus 33 
Three CSA-G 1296.2, 931.6, 846.2 31,24,13 0.58, 0.42 719 23.29 0.01655 0.9736 @ bus 18 
CSA-C 1061.3, 1215.1, 1047.2 14,2631 0.47, 0.53 99.6 29.83 0.00644 0.9823 @ bus 25 





4.2. Results Of Investigating THDv 
The aim of this section is to investigate the impact of adding the previous optimal calculated DG 
units to the RDS with the assumption that all added units are inverter-based PV units to get the worst 
condition of THDy. Detailed Simulink model of the studied RDS (IEEE-33 bus) is developed in phasor mode 
as it calculates the node and the branch voltage representing sinusoidal voltages and currents at a particular 
frequency. Firstly, the system is simulated without connecting any DG units, and without the presence of any 
harmonic-producing loads which ensures 0 % for THDy as expected at all buses. Then three other scenarios 
are investigated as illustrated in Figure 5. Such cases represent the three compromised solutions A, B and C 
using the modified CSA-G: 
a)  1* scenario: 3114.5 kW PV unit is connected at bus 7. As shown the voltage at all buses is distorted 
within the permissible limit. 
b) 2" scenario: 1344.7 kW PV unit is connected at bus 12 and 1164.25 kW PV unit is installed at bus 30; 
the results ensure that all buses are contaminated with harmonics within the permissible boundaries. 
c) 3 scenario: 1296.2 kW PV unit is added at bus 31, 931.6 kW PV unit is added at bus 24 while 846.2 
kW PV unit is connected to bus 13. As illustrated all buses voltage are distorted within accepted limit. 
As clearly shown, that the bus where an inverter-based PV system as a DG unit is installed has 
almost the highest THD, in most of the simulated cases. Besides, the distortion on other neighboring buses 
also increases, which means that neighboring buses are significantly affected; however the actual number of 
affected buses depends on the location of DG within the RDS. The THDy of the worst bus condition in all 
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studied scenarios results has increased to almost 1.97 % compared with the original case before DG addition, 
however it is still satisfying the limits of 3% according to the IEEE Std. 519™- 2014. Therefore, the 
achieved size and place of added DG units for the sake of improving the voltage profile and reducing the 
overall power system loss are accepted on the power quality issue of injected harmonics. 


3.0% 
---- One DG unit @ bus 7 


—— Two DG units @ bus 12 & 30 
— + =Three DG units at 31, 24 & 13 


2.5% - 


THDy (%) 


2.0% | 


1.5% 


1.0% |; 











Bus No. 


Figure 5. Impact of installing PV-based DGs (based on CSA-G) on THDy 


5. EXAMINING PROPOSED MODIFICATION FOR CSA 

To confirm the superiority of the proposed CSA-G and CSA-C algorithms over other algorithms in 
terms of quality of solution, the results of the original CSA algorithm and the modified one are compared 
against the results offered by other nature-inspired algorithms [35-36]. Such comparison is carried out to 
meet the goal of simultaneously minimize real power loss and improve voltage profile over IEEE 33 bus 
system using three DG units. The results are tabulated in Table 4, in which the performance of CSA, CSA-G 
and CSA-C and are compared with the documented results of PSO in [36], GA in [36], GA-PSO in [36] and 
TLBO in [35]. From Table 4, it may be noted that the proposed CSA-G offers a minimum value of power 
loss among all the algorithms referred there in. The voltage deviation obtained by the proposed CSA-G is not 
better than the original CSA but is better than the GA and PSO variants and the overall voltage profile of the 
system also gets improved, also CSA-G has the best economical index among all other algorithms. 


Table 4. Investigating Modified CSA Performance Over Other Methods 








; Optimal DG Loss reduction evaluation Voltage improvement 

Algorithm ; 

Ref. Size (MW) @ Bus Loss Reduction Ec. Voltage Lowest bus voltage 

(pu) (%) Index Index (pu) (pu @ bus) 

TLBO in [35] 1.1826, 1.1913, 1.1863 12,2830 0.1246 40.90% 41.20 0.005773 0.9829 @ bus 25 
GA in [36] 1.5, 0.4228, 1.0714 11,29,30 0.1063 49.62% 28.68 0.03511 0.9810 @ bus 25 
PSO in [36] 1.1768, 0.9816, 0.8297 8,13,32 0.1053 50.09% 28.26 0.03186 0.9806 @ bus 30 
GA-PSO in [36] 0.925, 0.863, 1.2 11,1632 0.1034 50.99% 27.66 0.01938 0.9808 @ bus 25 
CSA 1.251, 1.271, 0.611 14,24,30 0.087 58.76% 25.28 0.01834 0.961 @ bus 33 
CSA-G 0.846, 0.931, 1.296 13,24,31 0.079 62.55% 23.29 0.01655 0.9736 @ bus 18 
CSA-C 1.061, 1.215, 1.047 14,2631 0.0996 53.08% 28.83 0.00644 0.9823 @ bus 25 





6. CONCLUSION 

In this paper, a multi-objective optimization was applied to search for the optimal DG sizing and 
placing. DG is added in distribution networks for improving the voltage profile and reducing system loss 
keeping the injected harmonics values within accepted limits as a main power quality issue. The optimization 
problem is formulated from two conflicted goals, the total real power loss and the overall voltage deviation. 
A weighted sum method is presented to create the Pareto front and get the best compromised solution. The 
considered optimization problem is solved using modified CSA algorithm based on Gaussian distribution 
(CSA-G) for IEEE 33-bus radial distribution network in case of adding one or multiple DG units. Proposed 
economical and voltage indices are applied to check the achieved solutions of single, two or three DG units. 
A comparative study for CSA and the modified CSA-G and CSA-C for mathematical benchmark test 
functions depicts the better performance of the proposed CSA-G in terms of accuracy. It ensures that CSA-G 
is capable of enhancing the original CSA performance. Regarding the optimal DG sizing and location, it can 
be observed that the proposed method is able to find out the optimal location and size of DG, while, at the 
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same time, it reduces the real power loss and improves the voltage profile of the system. Upon achieving the 
compromised solution, the tested system is investigated for satisfying the voltage total harmonic distortion 
limits when incorporating PV based DG units according to the IEEE Std. 519™- 2014. Besides, the 
superiority of CSA-G performance is confirmed when compared against some other algorithms, hence, it can 
be deduced that the proposed CSA-G algorithm is a good choice over the original CSA and other reported 
algorithms for solving optimization problems. 
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